Understanding nanoscale temperature gradients in magnetic nanocontacts 
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We determine the temperature profile in magnetic nanocontacts submitted to the very large 
current densities that are commonly used for spin-torque oscillator behavior. Experimentally, the 
quadratic current-induced increase of the resistance through Joule heating is independent of the 
applied temperature from 6 K to 300 K. The modehng of the experimental rate of the current-induced 
nucleation of a vortex under the nanocontact, assuming a thermally-activated process, is consistent 
with a local temperature increase between 150 K and 220 K. Simulations of heat generation and 
diffusion for the actual tridimensional geometry were conducted. They indicate a temperature- 
independent efficiency of the heat sinking from the electrodes, combined with a localized heating 
source arising from a nanocontact resistance that is also essentially temperature-independent. For 
practical currents, we conclude that the local increase of temperature is typically 160 K and it 
extends 450 nm about the nanocontact. Our findings imply that taking into account the current- 
induced heating at the nanoscale is essential for the understanding of magnetization dynamics in 
nanocontact systems. 
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The discovery of the spin transfer torque P,i (STT) 
has provided the ability to manipulate the magnetization 
with an electrical current. This opened opportunities for 
new spintronic devices such as spin-torque operated mag- 
netic random access memories or nanosized spin-torque 
oscillators Q . Unfortunately, STT requires huge current 
densities, which can lead to substantial heating and early 
material fatigue. STT may also assist the magnetiza- 
tion switching in nanopillar geometries 0, H , and it can 
increase the magnetization thermal noise 16|, degrading 
the device performances. In nanopillar geometries, the 
magnetization dynamics takes place in a very confined 
region where the temperature is almost uniform Q , such 
that the temperature dynamics could be understood from 
simple experiments [5[. In contrast, the temperature rise 
and its spatial profile remain almost unexplored in 
nanocontact (NC) devices. This is problematic because 
magnetization dynamics takes place in a much wider area 
[9l.[lo| where the temperature might be very non uniform, 
potentially affecting the dynamics [H, 12|. In this con- 
text, it is essential to develop reliable tools to access to 
the local temperature during operation under large cur- 
rent densities. 

In this paper, we propose a methodology to access to 
the local temperature below the NC. The experimental 
bases of our method rely on the measurement and mod- 
eling of the electrical resistance of the nanocontact, com- 
bined with the study of the current-pulse-induced vortex 
nucleation [isj at different temperatures. The experi- 
mental results are compared with numerical results from 
tridimensional simulations of the heat generation and dif- 
fusion. The simulations allow us to access to the temper- 



ature profile. We conclude that at the very large current 
densities that arc commonly used 0, Qfor spin-torque 
oscillator behavior (around 4 10^ A/cm^), the tempera- 
ture increase is typically between 150 and 220 K with a 
very strong spatial gradient, which implies that it is es- 
sential to take into account the heating when aiming at 
understanding magnetization dynamics in nanocontact 
systems. 

We work on the NC depicted in Fig. [T] The top 
electrode forming the nanocontact is made of a Au(200 
nm)/Ti(10 nm) truncated cone with a radius at the base 
r„ = 60 nm imprinted in a 50 nm thick Si02 insulating 
layer (Fig [T}:) . The NC contacts the top of an extended 
spin valve (SV) of composition : seed layer 50/IrMn 
6/Co9oFeio 4.5/Cu 3.5/Ni8oFe2o 5/Pt 3 (thicknesses are 
in nm). The stack is grown on a 500 /^m thick substrate 
of intrinsic GaAs. The ends of each electrode (squares 
in Fig. [T^) are 250 /zm away from the NC, and they are 
electrically contacted with RF probes. 

The device magneto-resistance is typically 20 mO be- 
tween the parallel and antiparallel state for a total resis- 
tance in the parallel state of 6.2 17 at room temperature 
14| . Note that this value includes the resistance of the 
electrodes. 

Low temperature measurements (5 — 300 K) have been 
performed in vacuum inside a dark cryostat entirely sur- 
rounded by a radiative screen. The sample substrate is 
pressed on a copper thermal chuck, that is maintained at 
the temperature TappUed- The electrical probes and the 
radiative screen are also at TappUed- To measure the prob- 
ability of nucleating a vortex state using current pulses, 
we have used the protocol described in ref. |l3|: at a given 





FIG. 1. (Color online) (a) Top view of the sample. The 
bright areas are the gold electrodes. The substrate is cov- 
ered by 50 nm of Si02 and appears in black. The sample is 
contacted with electrical probes on the squares located at the 
outer edges of the image. At these positions and at the back 
side of the substrate the temperature is TappUed- The (col- 
ored) arrows show the direction of the heat flux away from 
the NC. (b) Scanning electron micrograph of the central part 
of the device, where we can see the 17 x 27 firn^ spin valve 
mesa. As the gold electrodes are deposited in a conformal 
way, the footprint of the buried NC appears as a circular dot. 
(c) Sketch of the NC. 
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FIG. 2. a) Square symbols: resistance of the parallel state for 
different applied temperatures. Line : linear fit. b) Evolution 
of the differential resistance with applied DC current for two 
different applied temperatures, 6 K (solid line) and 240 K 
(dashed line). The upper (dotted) curve is a parabola, used 
for comparison purpose only. The three curves are vertically 
offset for clarity. 
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applied temperature, we apply current pulses with am- 
plitudes 43 < / < 48 mA and durations Tp„;se=5-10 ns 
and measure the resulting microwave voltage spectrum 
to determine whether a vortex has been created [13|. The 
system is then reset to uniform magnetization, and the 
procedure is repeated 1000 times, with a waiting time of 
100 ms between each repetition, to get the vortex nucle- 
ation probability as a function of /, Tpuise and TappUed- 

Let us now summarize our main experimental findings. 

In a first step, we investigate the temperature depen- 
dence of the DC electrical properties (Fig. [2]). The zero 
bias resistance of the NC increases linearly with the tem- 
perature (^ = 6.4 mJl/K) from a residual resistance 
i?o of 4.3 f2. While applying some finite bias current, 
the differential resistance in the magnetic parallel state 
increases (fig. [Hd) almost with the square of the DC 
current for all investigated applied temperatures. This 
indicates that the dominant effect leading to the extra 
resistance 5R is Joule heating. The curvature of 5R is 
independent of the applied temperature. Surprisingly, 
this suggests that the current-induced temperature in- 
crease is independent of the resistance variations with 
the temperature in this system. Our understanding of 
this unexpected behavior is the following. Let us artifi- 
cially split the electrical system into two different parts: 
the Au electrodes on one side, and the nano-contact zone 
and the spin valve stack (NC-I-SV) on the other side. We 
consider that the Au parts are thick enough to be consid- 
ered as bulk materials, with negligible residual resistance 
^Octctztrodas ~ 0' ^hc total rcsldual resistance Rq is thus 
essentially that of the NC-f-SV part, written Rnc+sv- 



We believe that electron scattering due to alloy dis- 
order [isl in the spin-valve and interface reflexions 
are predominant in R^c+sv, which has two conse- 
quences, (i) The resistance below the NC is essen- 
tially independent of the local temperature (in practice, 
^Ij^pT and ^j^^^T « i?o)- (ii) The residual re- 
sistance i?o is mainly concentrated near the NC/SV con- 
tact. The latter point was confirmed [31 by studying the 
dependence of the total resistance at room temperature 
with the NC radius, which yielded an effective resistance 
area product {RA) of the NC of 20 mf^./xm^ at 300 K. 
This corresponds to an interface resistance Rnc = 2 
for the 60 nm radius NC system. As a consequence, the 
major part of the Joule losses occurs at the interface be- 
tween the SV and the NC and it yields a thermal power 
which is independent of the applied temperature. 

This should result in a temperature peak underneath 
the NC. The heat extraction is ensured by thermal con- 
ductivity through the gold top electrode, the rest of the 
SV stack in series with the gold ground electrode, and 
the GaAs substrate. Note that the thermal conductiv- 
ity of bulk gold is almost constant in the temperature 
window, such that we consider: '^""''^''"'^^ « 0. We 
also consider the thermal conductivity of the SV stack 
and substrate [itI as temperature-independent. Assum- 
ing that we can use an effective thermal conductivity K^f f 
relevant for our geometry and our material combination, 
the temperature increase below the NC {ST^c) is then 
independent of TappHed, and can be written as 
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TABLE I. Electrical and thermal conductivities used for 
the thermal simulations. The parameters of the substrate, 
Si02 insulating layer and SV are taken as temperature- 
independent. The SV is considered as a uniform material 
with average conductivities. 



where Ae/ / is an effective distance relevant for our geom- 
etry. Since the system presents a complex tridimensional 
geometry, simple considerations are irrelevant to deter- 
mine quantitatively Ae//. 

Instead, we now perform finite element simulations in 
order to evaluate the temperature rise and its spatial pro- 
file. Our solver (Comsol 4.2a) inputs the geometry of 
Fig. [T] and accounts for the electron transport [3| and 
the resulting Joule heating in the electrode, in the SV 
stack and in the NC with parameters listed in Table HI 
An intcrfacial electrical resistance is inserted between the 
NC and the SV, in line with our previous findings. The 
Joule-effect-related heat sources are the distributed re- 
sistive losses, plus this intcrfacial resistance. The heat 
diffuses in the entire volume, including the GaAs sub- 
strate and the insulating Si02 layer that surrounds the 
NC. 

The temperature profile in the stationary regime is cal- 
culated with the following boundary conditions : the gold 
leads terminate at the position of the contact probes, set 
at Tappiied- Thc substratc back side is also pinned at 
Tappiied- The remaining physical boundaries have free 
temperature conditions, with zero outgoing heat fiux. 

A first test of the simulation accuracy is to look at the 
temperature dependence of the resistance of the whole 
structure (Fig. The experimentally obtained linear 
increase of the resistance at warming (Fig. 2a) is reason- 
ably reproduced: the residual resistance i?o = 5.4 is 
1.1 f2 too large, and the simulated dR/dT = 4.3 mil/K 
slope underestimates the reality by 33%. Most of this 
discrepancy results from our neglecting the temperature 
dependence of the resistivity of the SV. 

A second test of the simulation accuracy is to look at 
the current-induced increase of the resistance. This is re- 
ported in Fig. for representative applied temperatures 
between 10 K and 296 K. The simulated R{I) curves 
are parabolic with a curvature almost independent of the 
applied temperature, recalling the experimental behavior 
(Fig. [^b), excepted near 10 K. We find also a parabohc 
increase of temperature due to Joule heating, almost in- 
dependent of the applied temperature, of the order of 170 
K for a DC applied current of 48 mA at the center of the 
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FIG. 3. (a) Resistance versus temperature (simulation result), 
(b) Symbols: (Color online) differential resistance versus DC 
current (I) for different applied temperatures between 10 K 
and 296 K (simulation result). Lines: parabolic fits. 
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FIG. 4. (Color online) Temperature profile, (a) in the plane 
of the free magnetic layer and (b) in the depth of the sample. 
Grey areas corresponds respectively to the region below the 
NC for (a) and of the NC for (b). Fig. a inset : temperature 
map in the (x,y) plane in a 800 nm by 400 nm area. The 
bottom left corner corresponds to the NC center. Fig. b inset 
: temperature map in the transverse plane (x,z) in a 800 nm 
by 400 nm area. The maximum temperature 360 K, is reached 
at the depth of the free magnetic layer below the NC center. 
The color contrast scale is from 280 K (blue) to 360 K (red) 
in both insets. 



NC/SV interface. 

Let us finally look at the simulated temperature pro- 
file, cither laterally in thc free magnetic layer or across 
thc SV thickness underneath the NC (Fig. E]). The tem- 
perature profile essentially preserves a cylindrical sym- 
metry around the NC (Fig. 2^ inset). It is peaked 5 nm 
below the NC (Fig. |4)d), i.e. right inside the free layer 
of the SV. The in-plane temperature distribution decays 
inversely with the distance outside of the NC (Fig. |3Ji) , 
with a width at half maximum A„,a,.,„ that varies from 
200 nm to 1200 nm when the NC radius increases from 
40 nm to 80 nm. If we take an average conductivity 
Keff = 15 W/(K.m), then the characteristic heat diffu- 
sion length Ae// is always greater than 1400 nm. Eq. [5] 
can then be viewed as a rule of thumb giving an upper 
bound for the temperature rise ST^^c in the NC. 

These conclusions are based on electrical measure- 
ments that do not inform directly on the temperature 
in the region of interest for magnetization dynamics. In 
order to double check our conclusion, wc shall now use 
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a different methodology to deduce the local temperature 
underneath the NC by the study of a thermally activated 
magnetization process occurring at that precise place, 
when a current is applied. 

Our probe of choice is the nucleation of a dynamical 
vortex state 0, [2^ [2lj as induced by pulsed currents. It 
has been shown previously [l3j that the vortex nucleation 
probability can be described by an Arrhenius law with 
a single activation energy Ea{I) which depends only on 
the total pulse amplitude. Here we shall measure the 
Arrhenius rate to get the real local temperature during 
the nucleation attempt. We use the classical description 
22 1 based on an Arrhenius- Neel law [i^, for which the 
mean nucleation time (Tnudeation) is : 



{^nucleation) — 'R) CXp 
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with T the local temperature and I/tq the attempt rate. 
The probability p of successful nucleation during the 
pulse (Tpuise) follows: 
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For each parameter set (/, Tpuise and TappHed), we would 
in principle need to determine the activation energy 
Ea{I), the attempt rate tq and the real local temper- 
ature T = Tappiied + STnc{I))- To reduce the number 
of free parameters, we make the following assumptions. 
The thermal equilibrium is reached in a time scale shorter 
than the pulse duration, the attempt rate is independent 
of J and T, and the temperature rise ((5T/vc(-f)) depends 



only on the pulse amplitude / and not on 



applied ; 



in line 



with our previous conclusions. We then use Eq. |l]to fit 
of the evolution of p with the applied temperature. 

Figure [5] represents the evolution of the nucleation 
probability with TappUed varying from 80 K to 296 K 
for 43 mA total current and a pulse duration of 10 ns. 

Two fitting procedures were implemented, and they 
are compared in figure [S] (i) In a first step we inves- 
tigate whether the data can be fitted disregarding the 
Joule heating (Fig. [5] dash (blue) curve ), i.e. consider- 
ing T = Tappiied- In that case the free parameters are 
Ea and tq. The fitting procedure used is based on the 
Levenberg-Marquardt method (23 | with no weighting of 
p. The results are summarized in Table HIl for the different 
parameter sets (Tpuise ,1) with the normalized residual of 
the fit. A mean activation time tq of 0.9 ns is found. Im- 
portantly, the heating-free fit cannot really reproduce the 
experimental behavior at low temperature, most proba- 
bly because the Joule heating induces a substantial error 
of absolute temperature in that range, (ii) In a second 
step, we account for the Joule heating. We consider the 
real temperature as T = TappHed + STNcil), and perform 
the fitting procedure with fixed tq at 250 ps and with 
only two free parameters Ea and STnc- This procedure 
is examplified in FiglS] (red line) . The results are listed 
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FIG. 5. Variation of — /n(l— p) with the applied temperature. 
Blue dashed line (color only): fit by an Arrhenius law without 
Joule heating. Red Line (color only): fit by an Arrhenius 
law taking in to account Joule heating and with a constant 
To = 250 ps. 
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TABLE II. Energy barriers Ea{I) and inverse attempt fre- 
quency To resulting from fits performed on the experimental 
data disregarding the Joule heating (i.e. T = TappHed)- 



in table mil The convergence of the fit is significantly im- 
proved for low temperatures and for high currents (i.e. 
low Ea{I)); i.e. in the cases where the role of the heating 
is the most important. 
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TABLE 111. Energy barriers Ea{I) and current- induced tem- 
perature rise resulting from fits performed on the experimen- 
tal data taking into account Joule heating with fixed to = 0.25 
ns. 

From these fits, the rise in temperature due to Joule 
heating can be written with a satisfactory agreement as 
STnc{I) = with r] 10^ K/A^. For instance, for 
a total current of 48 mA, the real temperature increase 
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during the pulse lies between 150 K and 220 K, which 
is consistent with the value of 170 K that wc predicted 
previously using thermal simulations. 

In summary, we have studied the temperature rise and 
its profile in magnetic nanocontacts submitted to the 
very large current densities that arc commonly used for 
spin-torque oscillator behavior. We accessed the tem- 
perature changes occurring during the application of the 
current using two experimental methods. They are: the 
modeling of the transport properties and, the modeling 
of the nucleation rate of vortex structures. Both methods 
were implemented in various current conditions and ap- 
plied temperature environments. We confronted our data 
to finite element simulations of the current and heat dis- 
tribution in a realistic tridimensional geometry. For the 
current densities typically used on nanocontacts in the os- 
cillator regime, the local increase of temperature reaches 
typically 150 to 220 K and it extends over an area with 
a typical diameter of 450 nm about a nanocontact of 60 
nm radius. 

In conclusion, thermal noise may be thus substantially 
higher than commonly expected, with detrimental con- 
sequences when using the nanocontact as spin torque os- 
cillators. In addition, this heating probably leads to sub- 
stantial changes in the magnetic properties, for instance 
in the exchange bias field of the reference layers. Finally, 
the temperature gradients arc unusually strong, and may 
lead to additional sources of spin torques [25[ that can 
play a role in the magnetization dynamics. Our findings 
imply that taking into account the current-induced heat- 
ing at the nanoscalc is essential for the understanding of 
magnetization dynamics in nanocontact systems. 
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